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This  paper  deals  primarily  with  3-D  effects  and  efficiency  enhancement  methods  in  a  steady  state 
FEL  amplifier  configuration.  We  treat  finite  transverse  dimension  effects  associated  with  i)  the 
wiggler  field,  ii)  electron  beam  and  iii)  radiation  beam.  Our  formulation  includes  efficiency  en¬ 
hancement  schemes  such  as  spatially  contouring  the  wiggler  field  as  well  as  accelerating  the  electron] 
beam.  Finally,  a  3-D  example  of  a  10.6  u  FEL  with  enhanced  efficiency  is  given. 
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THREE  DIMENSIONAL  NON-LINEAR  THEORY  OF  THE  FREE  ELECTRON  LASER 


Numerous  publications  have  treated  the  1-D  Free  Electron  Laser  (FEL) 
(1-10) 


As  of  the  writing  of  this  paper,  3-D  effects  in  the  FEL 

(11) 


mechanism. 

have  received  little  attention.  11 ^  It  is  the  purpose  of  this  paper  to  pre¬ 
sent  a  general  non-linear  3-D  formulation  of  the  steady  state  FEL  amplifier 
configuration  including  the  various  efficiency  improvement  schemes. 

When  a  cold,  axially  propagating  electron  beam  enters  a  physically  real 

( 12) 

w  ggler  field,  a  shear  in  the  axial  electron  velocity  results.  The  axial 

velocity  ..hear  across  the  electron  beam  is  equivalent  to  a  beam  temperature, 
which  if  large  would  significantly  reduce  the  fraction  of  electrons  trapped 
in  the  ponderomotive  buckets. 

The  radiation  beam,  which  experiences  both  diffraction  as  well  as  re¬ 
fraction,  can  be  represented  as  the  superposition  of  the  input  field  and  ex¬ 
cited  field.  In  the  absence  of  a  detailed  3-D  analysis,  it  could  be  argued 
that  these  two  fields  could  destructively  interfere  on  axis.  This  would  re¬ 
sult  in  a  decrease  in  the  depth  of  the  ponderomotive  potential  buckets  and 
could  cause  detrapping  of  the  electrons.  The  results  of  our  analysis,  how¬ 
ever,  indicate  that  destructive  interference  does  not  take  place.  A  number 

(7  9  po ) 

of  efficiency  enhancement  schemes  for  the  FEL  have  been  identified.  ’  * 

These  include:  i)  contouring,  spatially,  the  amplitude  and/or  wavelength  of 
the  magnetic  wiggler  field,  ii)  accelerating  the  electron  beam  by  applying  an 
external  D.C.  accelerating  electric  field.  Our  formulation  will  include  and 
show  the  equivalence  of  the  above  enhancement  schemes. 

A  generalised  pendulum  like  equation  for  the  phase  of  the 
particles  is  obtained.  The  wiggler  field  has  transverse  spatial  gradients 
as  well  as  an  amplitude  and  wavelength  which  are  arbitrary  functions  of  axial 

position.  The  infinitely  long  highly  relativistic  electron  beam  is  assumed  to 
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be  initially  cold  and  tenuous  enough  so  that  space  charge  effects  can  be  neg¬ 
lected.  The  inclusion  of  beam  temperature  and  space  charge  are  straightfor- 

(7  9) 

ward  and  have  been  included  in  1-D  formalisms.  ’  The  evolution 

of  the  radiation  field  governed  by  the  particle  dynamics  is  determined  as  a 

funct i on  of  x , y , z . 

Then,  simplifying  assumptions  are  made  which  result  in  an  analytic  expres¬ 
sion  for  the  field  on  axis.  Finally,  the  field  off  axis  is  numer¬ 
ically  evaluated  and  applied  to  specific  examples. 


The  generalized  linearly  polarized  wiggler  and  radiation  field  are  re¬ 
presented  by  the  following  vector  potentials 


"  A 

A  (y,z)  =  A  (z)  cosh  (k  (z)y)  eos(fk  (z')dz")e 
~w  w  w  °  w  y 

o 

uu  ^ 

A  (x,y,z,t)  =  A(x,y,z)  sin  (~z  -U)t+cp  (x,y,z))e 


where  A  (z)  and  k  (z)  =  2u  /l  (z)  are  the  slowly  varying  amplitude  and  wave- 
’  •  w  w  w 

number  of  the  wiggler  field,  &  is  the  wiggler  wavelength  and  A  and  op  are  the 

V 

slowly  varying  amplitude  and  phase  of  the  total  radiation  field.  We  also  in- 

A 

elude  an  external  accelerating  D.C.  electric  field  E  (z)  =  -  3d>  (z)/dz  e  . 

~ac  ac  z 

The  equation  (in  Lagrangian  independent  variables)  which  describes  the  rela¬ 
tive  phase  between  the  electrons  and  the  ponderomotive  wave 
can  be  shown  to  be  governed  by 


,  ,2  3  k  11/  d  $  it  2  /  r  . 

UL  =  lie.  +  _JL  +  JlEll  _  jejWc  d,  cosh  (k  ' 
2  2  dz  2  dz  ~2  2  t  3z  w  "  w' 


dz  '  dz 


Y  Ym  c 
7.  O 


2\  m  c 
'  o 


A  cosh  (k  y)cos  Pk  dz ')c 

XI  XI  w  u 


+  fk^A^  cosh  (k^y)  A(*x,y,z)  cos  If 


where  ji  =  ty  ( z  ,t  ,x^  ,y°)  =  J  (tu/c  +  k^(  z  ') )  dz u/t  +op(x,y,z)  is  the  phase. 


,  *!*■  V  *  f.V 


T7T-T 


T  =  t  +  J  dz'/v„,  y  .  y  V  ,  Y.  -  (1  -  TV)-1/2, 


o 


2  0-L 


Y  =  ( 1  +  (  e  A  (  z )  /  (m  c  )  )  )  ,  V  =  U)  /(uu/c  +  k  (  z  )  -  djf  /dz  +  d  j  /■.!:;)  /  ' 

ox  '  '~w  o  z  w  T 

,  z 

r*/  r*  ^  r>d 

electron  axial  velocity,  x  ~  x  +  3  k  sin  j  k  dz'  and  y  **  y  cos  K  z  ar 


OX  W 


w 


o 


the  zeroth  order  transverse  electron  coordinates,  3  =  v  /c , 

ox  ox 


p  2  -1 /? 

v  =  e  A  / (y  m  c)  is  the  wiggle  velocity,  y  =  (l  -  v  /c  ) 
ox  w  o  o  o  o 


is  the 


magnitude  of  the  total  particle  velocity  and  K  =3  k  /  >  2  .  The  Lagrangini: 

^  o  ox  w  ' 


independent  variables  t  ,x  ,v  are  the  time  and  transverse  ••fx.rdinat-s  of 

o  o  o 


tide  at  the  entrance  to  the  interaction  region,  i.e.  z  =  0.  Ir  c-t  tnininr 


we  have  used  the  fact  that  ui  »  (l  +  V  /c)  c^w’  x  comPone’''t'  of  ’--V  •- 
tron  momentum  is  (|e|/c)  (A  fy,z)  +  A„Cx,y,z,~) ) -e  ,  I A  ]  >  >  I  A_  I  and 

•  1  #v  X 


~  VD  ^  c*  Equation  (2)  clearly  shows  the  relationship  between  the  vario...* 


efficiency  enhancement  schemes.  In  the  equation  for  the  phase  if,  contouring 


j 

the  wiggler  wavelength  or  amplitude,  i,e.,  dk  /dz  or  dA  /dz  is  equivalent  tr 


introducing  an  accelerating  field,  i.e.  dty^/dz. 


The  axial  velocity  shear,  due  to  the  wiggler  gradient  is  giver:  by  AV 


shear 


=  c(3  k  y  /2)  ,  while  the  energy  shear  is  =  Y  (AV  ,  /c)m  c  .  The 

ox  w  o  shear  o  shear  o 


1/2  2 

initial  depth  of  the  trapping  potential  is  |e|$tr  =  2^  2Y0Y0Z®0x  moC 


In  order  to  trap  a  substantial  fraction  of  the  electrons  we  require 


AC  ,  <  |  o  I  <t>  ,  this  places  a  limit  on  the  electron  beam  radius  giver:  tv 

shear  1  'trap 


rb  < 


(Yk  )' 
6  w 


8'{  ?  y  1/2 


ox 


£r} 
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i/i. 


I  s) 


An  axially  symmetric  electron  beam  injected  into  the  wiggler  field  in  (la) 


experiences  a  periodic  focusing  in  the  y-direction.  The  focusing  wavelength 

,-l 


along  z  is  2n  K  »  i  .  A  comment  is  in  order  concerning  the  neglect  of 
o  w 

space  charge  waves.  It  can  be  shown  that  the  ponderomoti ve  term  in  (2)  (term 


proportional  to  A  A)  will  dominate  the  space  charge  term  if  the  beam  density 
w 


satisfies 


- -m 


9 


.  Ij  2  1 

n  <  <  ( k  "  Y  '  A  A )  ( 2  n  y  m  c  ) 

C  W  20  w  o  o 


2-2  2  2 

The  radiation  field  satisfies  the  wave  equation  (V  -c  d  /dt  )A 


n  e  *J  e  where  the  current  density  is  Riven  by 


'Mx,y,z,T. )  =  -  jj^.  J  dt  J*  dx^  J*  dyQ  9(xo,yo)6(x-x)6(y-7)6(t-t)  — 


A  (y,z) 


1  2 

where  or  =  (in  lel^n  /m  )  ,  n  is  the  density  on  axis  outside  of  the  inter- 

b  1  o  o  o 

act.  i . u:  rerion,  and  U ( x  ,y  )  is  a  function  which  describes  the  initial  electron 

o  o 

bear  profile.  The  integrations  in  (5)  are  over  all  initial  entrance  times  and 
l  !"e.; verse  coordinates.  The  radiation  field  in  (lb)  can  be  represented  in  the 
•i  •!•::.  «  i  x  ,y  ,t )  =  (2i)  ‘‘a(x,y,z)  exp  i  (aiz/c-OJt  )e  +  c.c.  where  a  =  A  exp(ip) 

is  Mu-  .•  implex  field  amplitude  which  is  a  slowly  varying  function  of  z.  Gub- 

r.t.  Lt.ut  inc  ( )  into  the  wave  equation  yields  an  equation  for  a(x,y,z)  which  can 
readily  t">  solved  using  Fourier  transform  techniques.  We  find  that  a  =  a.+a 


wn>*  re  a  , 


is  t. he  homogeneous  solution  (input  field)  given  by 


1  »  ■  —  2 
a  )  =  -  i  a(k  ,k  ,0)  exp  (i(k  x+k  y  -  -r — 7-  )  )  dk  dk  where 

1  .  x  y  £  x  y  2uu/c  x  y 

k‘  =  k‘  +  k'  and  a(k  ,k  ,0)  =  IT  a(x,y,0)  exp  (-  i(k  x  +  k  y))  dxdy.  At  z  =  0 

V  y  X  y  J  0  a,  X  y 

t  ).<•  t  t.ly  field  is  the  input  signal.  The  expression  for  can  be  shown  to 
rive  the  well  known  Gaussian  radiation  beam  modes.  The  particular  solution 
(excited  Hold )  is 

ju  "  z  2n/uudt®  ®  A' 

a.  ( x  ,.v  *z )  =  -  j ^  J'dxoi>dyo3(xo’yo)^7COSh  (kwy'} 

CO  O  -°°  -00  Y 

.  /  /  ~  ~  ^  v  2  4  uu/c 

l((x-x  }  +  (y-y  ]  )2T^PT  e-i(?'  -  cp')  (6) 


wh>T"  primes  on  quantities  denote  that  they  are  functions  of  z'.  The 


,*J  »AV>  V  ■  •• 


expression  for  the  excited  field  in  (6)  completes  the  formal  part  of  our 
analysis.  Equations  (?)  and  (6)  describe  self-consistently  the  non¬ 
linear  i-D  steady  state  EEL  amplifier  configuration. 


Considering  the  case  where  k^y^  <<  1,  x  and  y  can  be  replaced  by  x  and 
y  in  (6).  Assuming,  a  low  gain  situation,  i.e.  | a  j  >  |a  |,  and  a  plane 

O  1  c. 

wave  form  for  the  input  field,  we  can  take  the  phased  (z,tQ,xo,yo)  to  be 

very  nearly  only  a  function  of  z  and  t  .  Choosing  a  Gaussian  electron  beam 

2  2  2 

profile,  i.e.  Q(x  ,y  )  =  exp(-(x  +  y  )/r.  ),  (6)  reduces  to 
o  o  o  ob 


2,2  2  n/uu 

.  U),  /c  „ 

/  \  -ib  2  f 

ajr,z)  =  -  — -  J 


I  v  b  J  2rr/uu 
o  o 


dt  z  ,  „ .  /  z-z'  +  iz  \ 

o  p  ,  ,  .  ,  icp(r,z  f  o  \ 

3TJ  dz  A»u  )  e  l - — .,2  2) 

o  \  ( z-z  )  +  z  / 


.  U,  ,  ,  x  /Z'2'  +  izo  \  r2\ 

exp  -x  ^(z  ,to,-zo^-:::F77j7-J 


(7) 


where  2  =  r'  w/?c  is  the  effective  Rayleigh  length  associated  with  the  ex- 

o  b 

cited  radiation.  The  1-D  limit  of  (7)  is  obtained  by  letting  zq  or  r^ 
approach  We  will  limit  ourselves  at  this  point  to  a  constant  parameter 
wiggler  and  consider  only  an  external  accelerating  potential.  Furthermore, 
we  will  make  the  constant  phase  resonant  particle  approximation.  In  this 
approximation  all  particles  are  assumed  to  have  the  same  constant  phase,  1}  . 

n 

To  obtain  the  total  radiation  field  we  first  evaluate  a^(r,z)  under  the 
assumption  that  |qi  |  <<  1  (this  will  be  shown  to  be  valid).  Taking  a  plane 
wave  input  field  of  amplitude  A..^,  the  amplitude  and  phase  of  the  total  field, 
on  axis,  becomes 


A.(r  =  0,z) 


=  A. 


in 


2 

+  u.  A 
o  w 


tan  ^  (— )  cos 

z  R 

o 


(8a) 


(8b) 


»(  r  =  0,x)  =  -  a‘  (A  /A) 


“*  1  /  2.  \  .  'Y 

tan  ( — )  sm  * + 

"  n 


Z+Z 


4-)1/:  «*T. 


where  a  =  O).  r,  /2c  / y  ana  is  obtained  from  the  stationary  solution 
o  b  b  *  o  R 

of  (2). 


As  an  example  of  a  10.6  FEL  util izinp  a  nirh  paver  C0o  laser  beam  as 

an  input  field  we  choose  an  electron  beam  energy  of  2r  MeV  (y  -  80',  t^am 

current  of  I  =  h  A  and  beam  radius  (Gaussian  profile)  cf  r  =  0.2  rcr .  Ouch  a 

n 

beam  has  a  peak  density  on  axis  of  n  =  .1.1  x  lO^cm  ~  (jj.  =  2.0  x  1C '  "’sec  * ) . 

o  b 

The  const  ant  parameter  wi  cpler  is  taken  to  have  a  mapnitude  of  ip  =  f..O  KG 

and  wavelength  of  JL  =  C.1J  cm  which  pi  ves  A  -  f  . x  lo"1  statvclt.  The 
w  v 

_o 

wipple  velocity  is  y  =  2.6  x  10  c  an i  the  input  field  power  density  is 

taken  to  be  F.  =  l  x  lO^W/cirF  which  rives  A.  =  0.30  statvolt.  .  Note  that 
an  i  n 

the  inequalities  in  (3)  and  (8)  are  well  satisfied. 

Our  first  numerical  illustration  is  one  in  which  the  acceleratinp  poten¬ 
tial  is  zero,  the  stationary  phase  is,  therefore,  ^  =  -  n !'? .  The  particle 

enemy  remains -venstant  and  the  total  radiation  amplitude  and  phase  on  axis 
is  fiver,  by  (6a,  b).  To  obtain  the  radial  dependon  'e  of  the  radiation  amp¬ 
litude  and  phase,  (7)  is  solved  numerically  and  the  results  are  shown  in 
Firs .  ( 1 )  and  ( ;  ) . 


The  index  of  refraction,  in  this 

case,  is  rreat.er  than  unity,  n  =  1  +  ( c / uu  )uy>/3z  >  1.  The  input  field, 
therefore,  tends  to  focus  alonr  the  axis  and  <  ends  to  defocus  the  electron 
beam.  The  net  radiation  enerry  flux  alonr  the  z  axis  vinterrated  from  r  =  0 
to  r  =  »•)  is  constant  since  for  1 urpe  r  the  radiation  amplitude  is  less  than 
the  input  field  amplitude.  The  pain  in  the  radiation  amplitude  at  z  =  km 
is  maximum  on  axis  and  is  0.17.  The  maximum  value  of  7)  is  alonp  the  z  axis 
and  is  approximately  0.007  rni  which  certainly  satisfies  our  small  phase 
approximation. 


(> 


1 


• yv>  *-  • 


Our  next  illustration  is  one  involving  efficiency  enhancement.  An 


accelerating  potential  <(>^(2)  is  chosen  such  that  cos  =  0.3.  The  gain 

in  radiation  amplitude  on  axis  at  z  *  4m  is  0.185,  see  Fig.  1.  Since 

the  energy  gained  in  propagating  the  electron  beam  through  the  potential 

4>  is  converted  into  radiation,  the  efficiency  can  be  defined  as 
ac  m  .2 

j)  =  |  e  |  (<f  (z)  -  <f>  (0))/y  m  c2  =  -[  ^  )  ^  A  A<r  =  °*z')  cos  H' 

11  ac  ac  o  o  V  mocZ  /  o  °  w  R 

The  efficiency  at  the  end  of  z  =  4m  is  ~  3.6%.  Figure  (2)  shows  the 
phase  i f  as  a  function  of  z.  Notice  that  the  index  of  refraction 
n  =  1  +  (c/w)3s?/3z  for  large  z  is  less  than  unity  on  axis  (defocusinq 
of  radiation)  and  becomes  greater  than  unity  for  large  r  (focusing  of 
radiation) .  Equations  (8a, b)  are  in  excellent  agreement  with  the  above 
numerical  illustrations  for  r  =  0. 

The  authors  appreciate  useful  discussions  with  I.  B.  Bernstein  and 
W.  M.  Manheimer.  The  authors  would  also  like  to  acknowledge  support  for 
this  work  by  DARPA  under  Contract  No.  3817. 
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